Andreev-like reflections with cold atoms 
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We propose a setup in which Andreev-like reflections predicted for ID transport systems could be observed 
time-dependently using cold atoms in a ID optical lattice. Using time-dependent Density Matrix Renormalisa- 
tion Group methods we analyse the wavepacket dynamics as a density excitation propagates across a boundary 
in the interaction strength. These phenomena exhibit good correspondence with predictions from Luttinger 
liquid models and could be observed in current experiments in the context of the Bose-Hubbard model. 
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The rich physics described by Luttinger liquid (LL) theory 
im is normally associated with interacting one-dimensional 
(ID) electron systems such as carbon nanotubes or lithograph- 
ically defined quantum wires. However, exciting progress 
in cold atomic gases experiments O O has seen aspects of 
this physics realised in a new context IH [Sj. This not only 
promises observation of effects such as spin-charge separation 
in a clean system closely realising the theoretical mod- 
els, but also provides a new viewpoint on transport properties, 
which can be studied in the context of coherent wavepacket 
propagation. This connection is strengthened by the use of re- 
cently developed time-dependent density matrix renormalisa- 
tion group (t-DMRG) methods |6|, which allow the computa- 
tion of dynamics for physically realisable lattice models, and 
the identification of parameter ranges in which LL model pre- 
dictions can be observed in experiments. Here we investigate 
this analogy for systems described by the inhomogeneous LL 
model, in which the electron-electron interaction varies step- 
wise from essentially non-interacting to repulsively interact- 
ing Q. This model is used to describe the coupling of quan- 
tum wires to higher dimensional leads which act as weakly in- 
teracting electron reservoirs, and predicts rich boundary phe- 
nomena, including Andreev-like reflection, i.e., reflection of 
hole excitations. We show that analogous Andreev-like reflec- 
tions can exist for ID atomic gases in optical lattices, where 
the many-body dynamics are well described by Hubbard mod- 
els 1 8 1, and that these could be observed time-dependently in 
current experimental setups. 

Andreev-like reflections are predicted by an inhomoge- 
neous LL model with Hamiltonian Q (h = I) Hn = 
f {dx/2%) [u{x)g{x){Tinf + u{x){djc^f / g{x)], where <l> is the 
standard Bose field operator in bosonization and n its 
conjugate momentum density, [n(x),(|)(x')] = i?){x — x'). This 
model describes low energy excitations with speed of sound 
u{x). The parameter g{x) characterises interactions, with 
Q < g{x) < 1 for repulsive interactions, g{x) = 1 in the 
non-interacting case, and g{x) > 1 for attractive interactions. 
When a propagating electron-like charge excitation (or den- 
sity excitation) is incident on a boundary with g{x) = gi on 
the left of the interaction boundary and g{x) = gR on the right 




FIG. 1 : (a) A propagating excitation exhibits normal reflections (top) 
or Andreev Reflections (bottom) at an interaction boundary depend- 
ing on the relative interaction strengths on the two sides, (b) Obser- 
vation via bosons in an optical lattice in 3 steps: Preparation of the 
initial excitation using a superimposed trap (left); propagation of the 
excitation towards the interaction boundary, formed by coupling the 
atoms off-resonantly to an additional internal state (right); and de- 
tection via measurement of the atom density in a region between the 
location of the initial excitation and the interaction boundary. 



of the boundary, the strength of the reflections are quantified 
by a reflection coefficient y={gL-gR)/{gL^gR)' For y > 0, 
excitations are partly reflected and partly transmitted at the 
boundary (see Fig. [T^, top). However if y < then excita- 
tions are transmitted with a larger amplitude 1 — y which is 
compensated by the reflection of hole-like charge excitations 
with amplitude |y| (see Fig.[T^, bottom). This is analogous to 
Andreev reflection when an electron is incident on a normal 
metal-superconductor boundary: The electron forms a Cooper 
pair in the superconductor, and depending on whether its en- 
ergy is higher or lower than the superconducting gap, a par- 
tial or complete hole is reflected. Here there is no gap, and 
so the holes reflected are partial holes |7|. This phenomenon 
is manifest in several effects predicted for transport through 
quantum wires - such as oscillations of the nonlinear cur- 
rent voltage characteristics and the appearance of fractional 
charge excitations in the finite frequency current noise ||9l. 
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However, imperfections including contact resistance between 
quantum wires and the attached electron reservoirs have so 
far prevented these effects from being observed. In this sense, 
cold atoms in optical lattices would constitute an ideal physi- 
cal system in which Andreev-like reflections can be observed. 

We study dynamics on the lattice because the physics of 
atoms in optical lattices is well understood on a microscopic 
level fSl, and t-DMRG methods allow exact computation of 
the dynamics. We first investigate an extended Hubbard model 
with offsite interactions for spin-polarised fermions (or hard- 
core bosons), which corresponds in the continuum limit to a 
Luttinger liquid Q. The Hamiltonian is given by (/i = 1) 

H = -JY, cJcj^Y^^inifii^i +^8,-/T,-, (1) 

{ij) i i 

where Cf annihilates a fermion (or boson) on site /, / is the tun- 
nelling rate between neighbouring sites, nt = cjcf is the num- 
ber operator for particles on site /, Vi is the nearest neighbour 
interaction energy and £/ denotes the energy offset of site / due 
to external potentials. This Hamiltonian is valid for /, V^tt <C 
CO, with n the mean density, and CO the band separation. In 
the limit aVi/vp > 1, where vf is the Fermi velocity and a the 
lattice spacing, the connection between LL physics and this 
model is approximately given by gi = 1/ ^yl ^aVi/vp. Out- 
side of this limit, a connection can be drawn by fitting the 
Luttinger form for density-density correlations to numerically 
calculated results for the lattice model. 

Offsite interactions can be generated with Fermions, e.g., 
using polar molecules or by coupling to Rydberg states, or 
with hard-core Bosons by loading strongly interacting atoms 
into excited Bloch bands |[TOll . However, the natural exper- 
imental situation is to have short range contact interactions 
between atoms, as described by the Bose-Hubbard model in- 
cluding only on-site interactions, with Hamiltonian ih=\) 

^ = -^1^ ^Ih^llUMni - 1) (2) 

{iJ) ^ i 

Analogously to eq. ([l]), bt annihilates a boson on site /, / is 
the tunnelling rate between neighbouring sites, rii = b]bi is the 
number operator for particles on site /, Ui denotes the onsite 
interaction energy shift between two atoms and 8^ denotes the 
energy offset of site / due to external potentials. This Hamil- 
tonian is valid in the limit where /, Uifi <C CO. Note that in the 
limit \Ui/J\^\Ui/^i\ ^ 1 and with a mean filling factor /r <C 1, 
we can obtain the off- site interactions of the extended Hub- 
bard Hamiltonian, eq. ([T]) directly from onsite interactions. 
Restricting to the manifold of states containing only singly 
occupied sites, we obtain off-site interactions in perturbation 
theory as V^^eff = —J^/^i —fijUi^x . Below we also go beyond 
this limit in our numerical calculations. 

We begin by studying dynamics in the extended Hubbard 
model, eq. ([T]) before returning to the Bose-Hubbard model 
below. In analogy to the case in quantum wires, we would 
like to observe a density excitation propagating through the 
system towards a boundary in the interaction strength. This 



density excitation is formed by creating a local dip in the ex- 
ternal potential (e.g., using a focussed laser beam, see Fig.[TJ)), 

8,- = -8oexp[-(/-xo)V(2a^)] 0) 

Here, a denotes the desired Gaussian width of the local den- 
sity excitation, centred on an initial site, xq, and 8o will con- 
trol the depth of the potential and hence the maximum den- 
sity of the excitation. We choose the state for r < to be 
the ground state of the Hamiltonian ([T]), with 8o(^ < 0) > 0. 
At ^ = 0, we switch off the local dip in the external potential 
suddently, 8o(0 = 8o0(— 0^ leaving a Gaussian shaped den- 
sity excitation. The last term denotes a difference in potential 
between the left and right of the barrier, with the barrier func- 
tion F(x\ which beginning at site will be taken to vary lin- 
early across the width of the barrier, sites, as F(/) = 0, 
/ < xij\ F{i) = {i-xij)/Mij, Xb < i < Xb +M^; F{i) = 1, 
/ > Xh -\-Miy. The interaction Vi will vary in the same way as 
V; = Vl + {Vr — VL)F{i). Initially we will consider the bound- 
ary to be sharp, i.e., there is a step at a given lattice site x^, so 
that Vi = Vl for i < x^, and Vi = VrJ > x^. The parameter Er 
should be adjusted so that the density on each side of the bar- 
rier is approximately the same. Note that obtaining a constant 
initial density close to the barrier is only possible in the range 
IVl — as otherwise large oscillations are observed in 

the density near the boundary. 

The initial groundstate (with 8o > 0) and subsequent propa- 
gation are computed by quasi-exact imaginary and real time 
evolution respectively under the Hamiltonian ([T]). This is 
made possible by t-DMRG methods |6|, which are applica- 
ble to ID many-body systems where the Hilbert space can be 
expressed as the product of a chain of local Hilbert spaces. 
The state of the system is effectively written as a truncated 
Matrix Product state representation IH, in which % states are 
retained in each Schmidt decomposition of the system. In our 
calculations we performed convergence tests to ensure the ac- 
curacy of our results (see |11| for a general analysis of the 
accuracy of these methods), and estimate errors smaller than 
a few percent in the presented values. 

In Fig. [2^,b we show shaded plots of the density of atoms 
at different lattice sites as a function of time. From these plots 
we clearly see the propagation of the initial density excita- 
tion, which splits into a right-moving and left-moving excita- 
tion. The left-moving excitation is incident on the left-hand 
boundary of the complete system, and plays no role in the fol- 
lowing discussion. The right-moving excitation, however, is 
incident on the interaction boundary, resulting in reflected as 
well as transmitted excitations. In Fig. [2^, where Vl = and 
Vr = /, we see a normal reflection of the excitation, and both 
the transmitted and reflected excitations have a smaller den- 
sity than the incident excitation. In Fig.[2}3, we have Vr = —J, 
and we observe an Andreev reflection, where a hole excita- 
tion is reflected, corresponding to a lower density at the sites 
it occupies. Commensurately, a larger amplitude excitation is 
transmitted at the interaction boundary. 

We can quantify this process by defining the amplitude of 
a density excitation to be the total integrated change from the 
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background density over the sites containing the excitation. 
The reflection coefficient R naturaUy foUows as the ratio of 
amplitudes of the reflected and incident excitations. In prac- 
tice, we identify a group of sites between the original loca- 
tion of the density excitation and the interaction boundary as 
the measurement region. As shown in Fig. |2|:, the density in 
the measurement region increases and then decreases as the 
initial right-moving excitation passes, then we observe either 
an increase or decrease in the density resulting from normal 
or Andreev-like reflections respectively. If the measurement 
region contains some part of the initial excitation, then the 
background density is computed between the time windows 
where the initial and reflected excitations pass. The reflection 
coefficient can be found as the ratio of the peak values in each 
of these time periods. Note that this definition can be used 
operationally in an experiment, where the integrated density 
over several sites can be measured, e.g., using flouresence or 
phase-contrast imaging with a focussed laser (see Fig.[TJ)). 

In Fig.[2]l we plot the reflection coefficient R for —J<Vr< 
/, choosing Vl = 0, 8o = 2/, and a = 3, and see clearly the 
crossover from Andreev to normal reflection. Note that the 
quantitative values (though not the behaviour with varying Vr) 
are dependent on the size of the initial excitation, as shown 
in the inset of Fig. |2|l. In order to compare our results with 
the known analytical result from LL theory, we extracted ap- 
proximate Luttinger parameters geff for the ground states of 
our extended Hubbard model, eq. ([T]), at the background den- 
sity value for each Vr. This was done by computing density- 
density correlations, and fitting the standard form for these 
correlations in a Luttinger liqud 1 1 1, 



{nofir) ^ -7^-2:^2 +^cos(2mr) ( - 



(4) 



where A is a constant, and n is the mean occupation per lat- 
tice site. Using these values we computed approximate re- 
flection coefficients, which are plotted in Fig. [2]l. They show 
very good quantitative agreement with R in our simulations for 
Vr > 0, and only small deviations for Vr < 0, especially for 
initial density excitations containing approximately one par- 
ticle. This agreement is better than might be expected, and 
demonstrates both the generality of LL theory in its applica- 
bility to low energy states of ID systems, and its continued 
applicability when we introduce by hand the boundary in the 
interaction strength. For comparison we have also computed 
geff based on the approximate analytical expression. 

Until now we have considered only sharp barriers with an 
immediate transition from Vl to Vr. We now investigate the 
time propagation of the excitations when Vi varies linearly 
over Mh sites, which provides a more realistic treatment of 
barriers that might be created in a real experiment. An ex- 
ample shaded plot of the density of atoms at different lattice 
sites as a function of time for an excitation exhibiting Andreev 
reflection with Vl = O^Vr = —J is shown in Fig.|3^. We see 
clearly that the extended length barrier spreads the resulting 
reflected and transmitted excitations in space (cf. Fig. 
However, the total amplitude of the reflected wavepacket is 
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FIG. 2: Numerical simulation results for propagation of an initial 
density excitation (80 = 2J, a = 3, xq = 45) across an interaction 
boundary from = to varying Vr at Xfy = 90 (Mfy = 1), with open 
boundary conditions, (a) Shaded density plot for Vr = /, showing a 
normal (positive density) reflection at the boundary, (b) Shaded den- 
sity plot for Vr = —7, showing an Andreev (negative-density, or hole) 
reflection at the boundary, (c) Difference from the initial value of the 
integrated density over sites 60-75, showing initial peak due to prop- 
agation of the initial density excitation, and a secondary maximum 
(Vr = 7, dashed line) or minimum (Vr = —J, solid line) due to the 
reflected excitation, (d) Reflection coefficients as a function of Vr, 
showing the comparison between simulation values (solid line), and 
estimated parameters from LL theory using numerically computed 
Luttinger para meters (from eq. (|4|, dashed line) and the analytical 
form g = 1/ y^l-\-aVi/vF (dotted line). The inset shows the depen- 
dence on the depth 80 with a = 2 (solid lines) and a = 4 (dashed 
lines) for Vr = —0.5/ (upper curves) and Vr = —J (lower curves), 
computed from measurement sites 50-70. 



actually increased in this case. This can be seen from Fig.jSJ), 
where we plot the reflection coefficient 7? as a function of 
Mh for different amplitudes of the incident excitation. This 
brings the reflection coefficients closer to the values obtained 
from the effective LL parameters. The finite- width barrier also 
smoothes the local density minima or maxima that appear at 
the boundary in the ground state (see Fig.|4^). 

To provide a simpler experimental implementation, we con- 
sider the dynamics of excitations at an interaction boundary in 
the Bose-Hubbard model. In the perturbation theory limit, it 
follows from \J/U\ <C 1, that |V,,eff| < /, and thus the am- 
plitude of all reflections will be extremely small. However, 
using numerical simulations we can treat this system exactly 
also beyond the limit in which perturbation theory is valid. In 
Fig. we show results from numerical simulations in which 
a density excitation (created as for the extended Hubbard 
model) propagates across a interaction boundary at with 
Ui<xh = 10/ and Ui>x^ = /. We see very clearly that a hole 
excitation is reflected, and reflection coefficients as a function 
of y^^eff = 1 /Ui>xb are plotted in Fig.|4t)., and compared with 



results for the extended Hubbard model with V, 
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FIG. 3: Reflections of a density excitation (80 = 27, xq = 45) from 
a thick boundary with the interaction strength varying from Vi = 
to Vr = —J Unearly over lattice sites, (a) Shaded density plot 
showing reflection of an initial density excitation with a = 3 from a 
boundary with = 10, showing the considerably broader reflected 
wave produced by the thick boundary, (b) Reflection coefficient as a 
function of barrier thickness (measurement sites 55-85), showing 
an increase in the amplitude of the negative density reflection as the 
barrier is increased, for initial wavepackets with a = 3 (solid line) 
and a = 2 (dashed line). The amplitude is also slightly larger for the 
narrower initial wavepacket. 80 = 0) = 2J. 




FIG. 4: Realisation of Andreev reflections within the Bose-Hubbard 
Model: (a) The density (fii) at each lattice site at time tJ = 5 (dashed 
line) showing the two counterpropagating density excitations, at time 
and ^7 = 30 (solid line), showing the transmitted and reflected excita- 
tion with Ur = J. Note the clear negative-density excitation centred 
near site 70. (d) Reflection coefficients for varying Ur, plotted as a 
function of Vr^q^ for the Bose-Hubbard model (solid line). These are 
compared with results from the extended Hubbard model (dashed 
line) and estimated parameters from LL theory (dotted line) com- 
puted as in Fig. 2. Parameters used were 80 = 27, a = 3, = 45, 
Mjy = I, Ul = lOJ, measurement sites 55-85. 



Effects that go beyond the validity of perturbation theory ac- 
tually yield a slight increase in the amplitude of Andreev-like 
reflections, and bring these results close to those predicted us- 
ing the estimated Luttinger parameters that were plotted in 
Fig. [2]l. In addition, we can use values of U ^ J, where 
IVeff/-^! > 1 without adverse boundary effects that prevent us 
from obtaining a smooth background density near the interac- 
tion boundary. Thus, it is possible to observe even larger am- 
plitude Andreev-like reflections in the Bose-Hubbard model 
than in the extended Hubbard model. 

In an experiment, the change in onsite interaction strength 
could be engineered simply in several ways. For example, 
lasers focussed on one side of the system could couple the 
atoms off-resonantly from their internal state, |a) to an addi- 
tional internal state |P). If a Feshbach resonance |12 | exists 
between two atoms both in state a so that their interaction 
Uaa is much larger than that between atoms in the two differ- 
ent states [/otp, then adding an admixture of the state P will 
reduce the onsite interaction strength in the region where the 
internal state of the atoms is |\|/) = |a) +(22 |P), where a\ 
and a2 are complex coefficients. Such a laser coupling could 
reasonably be focussed so that the coupling varies on a length 
scale of < 5^m, or approximately 10 lattice sites. Similarly, 
the initial density excitation could be prepared using a laser 
focussed over ~ 10 lattice sites. 

We have shown how time-dependent wavepacket dynamics 
of Andreev reflections in Hubbard models closely match the 
behaviour expected from LL physics, and how these reflec- 
tions could be observed time-dependently with cold atoms in 
optical lattices. This work could be extended to multi-species 
models, such as fermionic Hubbard models, and to other ge- 
ometries such as Y-junctions |13|. The boundary physics 
in each case has different characteristics which could be ex- 
plored time-dependently in the experiments. 
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